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Bistatic vs. Monostatic Measurement 


¢ Bistatic — Transmitter and Receiver * Monostatic — Transmitter and 
are at separate locations. Receiver are co-located. 

¢ Tx and Rx need a common time / ¢ Measures the Critical Frequency 
frequency reference. (vertical incidence) above the 

* Measures MUF of the path sending/receiving location. 
centerpoint between Tx and Rx. ¢ Tx & Rx can share reference. 

© Faye higher than F yitical- ¢ Transmit signal ‘buries’ the receive 
* Modest Rx dynamic range needed. echo. 
Due to T/R separation ¢ Large Rx dynamic range needed. 
lonosphere isaaepiel 


F-layer 


Fay 


E-layer 


Transmit ma ie Time / Freq. Reference 9 Receive Antenna 


Transmit Dipole Receive Small Loop 


Earth Earth 


Monostatic System Concerns 


¢ Operating full-duplex — reflection delay is about 1.6 
milliseconds. 
¢ Transmit signal is much stronger than receive echo. 
¢ Rx overload, Tx & Rx Phase Noise, Rx ADC 
performance are all important. 
¢ Hermes: works well. 
¢ Need high relay isolation for full duplex echos. 
¢ Leakage: Separate Rx Antenna relay > T/R relay is an issue. 
¢ Need to know precise Transmit time. 
¢ Equipment delay uncertainty impacts measurement. 
¢ Ethernet Packetization, Operating System scheduler, etc. 
¢ Need good frequency control. 
¢ Doppler shift is tenths of a Hertz. 


lonosphere 
© Fojitica, : Maximum frequency that is reflected at vertical incidence. 


¢ Faye : Maximum Usable Frequency that is reflected at angle 9. 
*  Faur = Foritica / Sin 8 


¢ Height of the F-layer > Reflection delay time. 
* Doppler Shift : Due to vertical velocity of Five 
* Ejayer : attenuates transiting signals (twice). 


F critical Fur 
(8 = 90 degrees) (8 < 90 degrees) 


At lower incidence angle, the wave 
travels a longer interaction path 
through the ionosphere. 
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Linear Chirp Approach auiiill il 
\ al Hl I ill ih Ill 
¢ Transmit a linearly swept FMCW waveform. —-—-. = ; 
¢ Measure instantaneous frequency offset or time offset 
between transmitted and received chirps. 


¢ While transmitting (i.e. Full Duplex). 


lr 


Time delay due to 
ionospheric height, 
no Doppler shift. 

Infer the Range or 
Height of the layer. 


frequency 


Frequency delta 

due to ionospheric 
height, no Doppler 
shift. 


time 


Doppler Shift > Range Error 


¢ Vertical velocity of the F-layer causes Doppler Shift of Received chirp. 


¢ Doppler shift causes range error in the measurement. 
¢ Range error is equal and opposite between up-chirp and down-chirp. 


(+) Doppler shift due to 
falling ionospheric layer 
(approaching earth) 


Causes (-) Range 
error due to (+) 
Doppler shift 


Causes (+) Range 
error due to (+) 
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Digital Signal Processing 


¢ Difficult to see a received echo buried deeply in transmit signal. 


¢ Approach: Correlation of received signal against transmitted 
signal replica — “Matched Filter”. 


¢ Strong peak due to transmit signal — helps identify precise transmit time 
reference. 


¢ Weak peak due to received echo — compute contribution of received 
signal across the entire chirp time. 


¢ Lots of taps > good dynamic range, good time resolution. 
¢ Multiple sweep integration — reduces noise as sart(N). 


¢ Implementation: FIR filter with time-reversed taps. 
¢ Taps: Time-reversed transmit Sequence. 
¢ Problem: 10° FIR filter taps > 10'* computations. 
¢ Solution: Frequency domain filter. 10° filter taps > ~10’ computations. 


¢ FFT (both signal and taps) > Filter (multiply Rx’d signal by Filter taps) 
> IFFT 


oystem Block Diagram 
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Time Domain Windowing 
¢ Discontinuities at the edge > spectral leakage. 
¢ Obscures the echo. 


¢ Solution: ‘window’ the received sequence in time. 
¢ Blackman-Harris window. 


¢ Implement as lowpass filter preceding correlator. 


¢« Extreme Positive and Negative frequencies are reduced in 
amplitude. 


Low-pass filter amplitude response 


amplitude (4) amplitude / Xx 
frequency = frequency 


-) = frequency below channel center 
(+) = frequency above channel center 


Windowing Reduces Correlation Sidelobes 


¢ Transmit peak near zero milliseconds, vertical axis dB. 
Blue — no LPF Red — Blackman-Harris LPF 


Top Block 
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Time (ms) 
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Gnuradio Flowgraph 


Options 
ID: top_block 
Generate Options: QT GUI 


ChirpSource 
Chirp Rate: 15 
Sample Rate: 48k 
Direction: 0 


Amplitude: 950m 
Deviation: 20k 

Zero Pad: 1 

Window duration: 3m 


Complex to Mag*2 


Variable 
ID: samp_rate 


Value: 384k 


QT GUI Chooser 
ID: Frequency 
Num Options: 2 
Default Value: 3.7M 
Option 0: 3.7M 
Label 0: 80m 
Option 1: 7.025M 
Label 1: 40m 


hermesNB 
Revr 0 Frequency, Hz.: 3.7M 
Revr 1 Frequency, Hz.: 7.2M 
Transmit Frequency, Hz.: 3.7M 
Rx Sample Rate: 384k 
Rx Preamp Off/On: 1 
PTT On Mutes Rx: 0 
PTT Off Mutes Tx: 0 
Tx PTT mode Off/Vox/On: 2 
Tx Drive Level (0..255): 25 
Ethernet Interface: eth0 
HPSDR Clock Source: OxF8 
Alex Tx Antenna: Tx1 


Alex Rx Antenna: Tx Ant via T/R Relay 


Alex Rx HPF: AutoTrack RxOFreg 
Alex Tx LPF: AutoTrack TxFreq 
Verbose (1=on, 0=off): Off 
MAC Address or *: * 


File Sink 
File: out_filename 
Unbuffered: Off 


Variable 
ID: prefix 
Value: /home/tom/Desktop/ 


QT GUI Entry 
ID: TxDrive 


Import QT GUI Entry 
Import: numpy ID: Alpha 
Default Value: 125m 
Import 
Import: datetime 


Default Value: 25 


Variable 
ID: out_filename 
Value: /home/to...22.16.15.32 


Windowing LPF 


Low Pass Filter 
Decimation: 1 
Gain: 5 

Sample Rate: 384k 
Cutoff Freq: 10k 
Transition Width: 7.5k 
Window: Blackman 
Beta: 6.76 


Correlator 


FFT Filter 
Decimation: 1 
Taps: numpy.fromfile(‘/hom... 
Num. Threads: 1 


Append file: Overwrite 


Post Processing Software 


¢ Python program 
¢ Reads separate received up-chirp and down-chirp files. 
¢ Isolates each individual sweep. 
¢ Time-aligns transmit peaks to zero time. 
¢ Discards ‘dud’ sweeps (sweep #0, corrupted, etc.) 
¢ Performs non-coherent integration. 
¢ Separate integrations for up-chirp and down-chirp 


¢ Plots output 
¢ Up and Down chirps with different colors. 


Spyder? : a good Python interactive development tool. 
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Results 


¢ Excellent performance of Hermes module in full duplex 
mode — phase noise and ADC performance yield about 
90~100 dB of post-correlator dynamic range. 


¢ F-layer measurements done on 80m. 
¢ Daytime sweeps yield null results due to excessive E-layer 
attenuation. 
¢ After the sun sets, the E-layer de-ionizes. 


¢ F-layer Echos clearly seen. 
¢ Occasional double-transit echos seen. 
¢ Occasional O-wave and X-wave are resolved. 


3.6 MHz — Local Noon 


¢ No F-layer reflection seen. E-layer attenuation high at 
3.6 MHz. 


¢ Spurious responses at 3.07 msec and 4.00 msec. 


80M Noon F-layer 
‘en Transmit signal correlation peak set to 0.0 msec. 


— Up Chirp 
— Down Chirp 
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3.6 MHz. Evening 


F layer echos near 1.7 milliseconds. -0.38 Hz Doppler shift. 
Double-transmit echos near 3.4 milliseconds 


Height: 254 km. Vertical Velocity: +15.4 m/sec. Artifact at 4 milliseconds 


120 
— Up Chirp 
> | — Down Chirp 


F-layer echos near 1.7 milliseconds 


magnitude, db 
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Ordinary and Extra-ordinary Waves 


¢ Alinearly polarized signal decomposes into RHC and LHC 
components in circularly birefringent ionospheric layer. 


¢ Earth’s DC magnetic field biases electrons. Causes electron spin 
precession (Electron gyro-resonance about 1.4~1.8 MHz). 
¢ Called O-wave and X-wave. 


¢ O- and X- Independently reflected by F-layer. 


aaa 80M Evening F-layer 


— Up Chirp 
— Down Chirp 


Ordinary: (O) & Extra-ordinary (X) F-layer echos 
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Further Work 


¢ E-layer measurements require 160m transmit antenna 
of reasonable efficiency and match. 
¢ With good zenith radiation pattern. 
¢ Identifying the O-wave and X-wave requires two 
receive antennas, two coherent receivers, etc. 


¢ Baseband processing can synthesize RHC and LHC froma 
pair of linearly polarized receive antennas. 


¢ Thanks to: 


Andrew Martin VK30E, John Petrich W7FU, and Phil Harman 
VK6PH for review and comments. 
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